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Velocity relative to a blade row 
Induced velocity defined by eqn. (2.37) 
Variable, equal to — -v-- 

f iC 

f Variable, defined by eqn. (2.l6) 

G Variable, defined by eqn.(2. 1l) 

H Total enthalpy 

I,J lUnctions defined by eqns. (2.18) to (2 

i,3 Index number for grid points 

£ Total length of camber line 

m Summation Index 

N Number of blades 

Q Volume flow rate 

r Radial coordinate 

S Length along any coordinate direction 

T Rinct ion defined by eqn. (2.26) 

t Half the thickness of blade 

u Absolute velocity in a-direction 

V Absolute velocity in 3-direction 

w Absolute velocity in e-direction 

X Variable, equai.-t© N'T 

z ^ 

Greek Symbols 

3 Meridional streamlin^^^ coordinate 
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r Circulat ion. 

Y Vortex strength 

'5 Dirac’s delta function. 

^ tan (slope of the camber line) 

? Vorticity 

n Coordinate normal to the camber line 

0 Angle, defined by eqn. (2.22) 

0 Circumferential coordinate 

Correction factor for v^, defined by eqn.(2.50) 
Parameter equal to ml 
C Coordinate along the camber line 

^ Source strength 

T iUnction of 0, defined by eqn. (2.9) 

$ velocity Potential 

"P Stream function 

CO Angular velocity of the rotor 

Subscripts 

0 At the singularity 

1 Upstream of singularity 

2 Downstream of singularity 

+ Suction side 

_ pressure side 

* ' Por flow with infinity large number of blades 

h At hub of the turbomachine 



i Iraet edge 

jc Blade surface 

H Trailing edge of the blade section 

m Value corre spend jrg to summation index m 

p* Axisymmetric potential flow 

q Value related to sotH-ces/sinks 

t At casing of the turbomachine 

a,e,e In the direction of respective coordinates 

Y Value related to vortices 

n In the direction normal to the camber line 

5 Along the camber line 

Superscripts 

“ Pertairbation value 

- Value relative to a blade row 
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Derivative w.r.t. 3 
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SYUOPSIS 

A MJMERICAL SOIUTIOKf OF THE IITVERSE PROBIEM FOR HYDRiOLIC TaRBOMAGHElES 

US mo A SINGULARITY METHOD 

N.K. JAIN 

MASTER OR TECHNOIDGY 
Department of Mechanical Engineering 
Indian Institute of Technology, Kanpur. 

In the present: work a computer program and computational 
results are presented for the blade design of a hydraulic turbomachine, 
using a singularity method based on Kashiwabara’ s work. Blade shapes 
are obtained £rom the specified velocity distributions on blade surfaces 
and camber line lengths on each stream surface. The flow thrcu^ the 
turbomachine is assumed to be incompressible and invisoid. 

Fortran program listing and description of the computer 


programs are also given. 



CHAPTER I 


IHTRODUCT lOH 


A designer is basically concerned with two essential aerodynamic 
problems, Mrst is the "Analysis problem" in vdiich he is concerned 
with the accvirate prediction of the design and off-design performance 
of a given ^st®L of blades. Ihiie in the "Inverse problem" the 
designer requires to know what cascade of blades will perfoim a 
prescribed duty. In hydraulic tiirbomachines for given hub and casing 
geoaetry, he is at liberty to choose some among the following* 
tiie position of inlet edge of the blade* camber line lengths at 
various stream siirfaces of revolution; camber line shapes; velocity 
or pressure distributions on sxiction and pressure surfaces of blade 
sections on a number of stream sviifaces; thickness distribution; 
depending upon the method used for the solution. 

A theoretical approach to the aerodynamic design of annulus 
and blade cascade of turbomachinery presents a formidable problem. 

It involves the solution of partial differential equations of the 
fluid flow in complicated geometries. It is only with the help of 
large fast computers that a realistic attempt can be made at all. 

In its entirety suda an attempt woiild involve the solxition of the 
eqimtions of the three-dimensional, unsteady, ccsnpressibie, viscous 
flow throu^ the complex of rotating and stationary blades and asso- 
ciated ducts of a turbcmachine. She size and speed of the present 
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day ccmpiiters is not adequate for this purpose. For is the state 
of o-ur knowledge about turbulent, separated and secondary flows. 

Usually efforts ha-ve been concentrated on the formulation of meaningful 
combinations of two-dimensional models of the processes involved to 
model the real three-dimensional flow. All references to three- 
dimensional methods should hereafter be understood to mean such 
formulations . Fmerical techniq'ues are applied to solve the 
associated two-variable non-linear partial differential equations 
which result from these models. The ultimate goal is a series of 
computer programs throu^ v\Jaich an aerodynamic design can be formulated 
and its off-design performance accurately predicted. 

The present numerical solution is a computer program for the 
design of blades of ajqy single stage hydraulic turbcmachine. The 
computer may be regarded as a simulator and is used in the sane manner 
as an experimental rig of great flexibility in which the effect of 
changes in design parameters may be investigated rapidly. A design 
can be finalised with some ceidiainty before actual manufacture and 
testing with their problems of hi^ cost and time delays. 

While designing the blades of any axial, radial or mixed 
flow turbomachines, it would be desirable to obtain blade shapes 
by prescribing velocity distributions along the blade surfaces 
as to avoid botmdary layer separation, cavitation and other problems. 

liaa^ fo^ the inverse method for two-dimensional^^^^^^^^^^^^^^^^^^^^^^^^^^ 

and three-dime nsional cascades are available, Stanitz^ gives a 
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two— dimensional approximate method widely -used for getting an 

approximate design of blade cascade exclijding tte nose and tail 

2 

region of blade element profile, Wu published a method in which 

and velocity distri 

the shape of the streamline between blade element profiles^ along 
this stream line are prescribed together with the blade spacing, and 
the resultant blade element profile and velocity distribxiiion along 
it are determiined. The method involves approximations that are 
accurate for hi^ solidity cascades, but has no direct control over 
the velocity distribution along the profile surface. St anitz and 

3 

Sheldrake discussed a two-dimensional design method pa?escrlbing 

surface velocity distributions as a function of arc length of camber 

line for hi^ solidity cascades, but it can not determine the blade 

shape in the region of stagnation points and also it is not applicable 

to blade design of turbcmachines in which meridional stream lines are 

4 

curved, Wilkinson developed an improved version of singularity 
methods developed by Martensen and by Jacob and Riegeis. It needs 
initial blade configviration, chord lengths, mean flow angle and 
thickness distribution in addition to the pressure or velocity distri- 
butions on the suction surfaces. Moreover, velocity distributions 
on pressure suif aces can not be specified in this method. The method 
works only for the design of blades of large pitch/chord ratio and 
stagger angles of less than 60°, In Eef. (g) Y/ilkinson has given 
some guide lines for a three-dimensional design method, vhich is 
similar to the above and is a simple modification of the 
method described there, Murugesan has given an improved version 


4 

of Wilkinson's method, which operates by the systematic modification 

of geometiy of an arbitrary initial blade shape to satisfy the 

reqxhrement of desired velocity distribution, but is limited to two- 

- rr 

dimensional incompressible flow. Reilly ' has given a quasi-three- 

dimensional singularity method which approximates the stream surfaces 

by a connected series of truncated cones. Mural' s method is a 

three-dimensional method in which circulation distribilion is prescribed 

9 

It is restricted to thin blades. Kashiwabara has given a three- 
dimensional singularity method for the design of blades of hydraulic 

a 

turbomachines, which is an improvement over Mural's method , This 
method appears to be very useful for hydraulic turbomachines of almost 
any kind, and thus is selected for the numerical solution presented 
here. 

At this point, the basic question one can ask while starting 
the solution for an "Inverse problem" is how to specify an optimum 
velocity distribution on the sections of blade element surfaces. 

Of course, the answer is not very easy, but a number of investigations 
(16, 12, IS, 14 ) have been made which attempt to give guide lines 

for the optimal velocity distributions for two— dimensional subsonic 
cascade problems. The optimization is with respect to the following 
considerations. 

i. Minimization of boundary layer separation losses. These 

pan be avoided by prescribed velocities that do not decelerate 
rapidly enou^ to cause separation. 
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2. Avoidance of cavitation in incaapressible flew, and shock 
losses in compressible flow. Ere&cribed velocities should not 
exceed certain maximua values dictated by these pheneneaa. 

5. ihr compressible fLov/ the desired flow rate can be assured by 
prescribed velocities that do not result in "choked flow" conditions. 



CHAEPER II 


THE AIAIYTICAD IDMUIiATIOH OH THE EROHLiM 


lii the present solution of the inverse problem, the flow 
throu^ the turbanachine is assuaed to be incompressible 
and inviscid. Further, since the flow far upstream is irrotational, 
it will remain irrotational everywhere with respect to an inertial 
frame (lab, fraae). 

First the meridional potential duct flow solution is found i,e. the 
flow without blades or blade forces. Stream surfaces, even with 
blades present, are assuaed to be close to the stream surfaces of the 
p)Otential duct flow without blades, defined by the hub— casing boundaries. 

Blade shapes are found on these potential flow stream surfaces and 
the flow on such surfaces, even with blades present, is assumed to be 
a two-dimensional flow. Yortices, sources and sinks are used to replace 
the blade in these two-dimensional sheets, analogous to plane potential 
flow singularity methods for solving the inverse problem. 

When the blades are present, the velocity potential inside the 
blade row is assimed to consist of two parts; one due to an infinitely 
large nmber of blades, and the other part, a perfurbation velocity 
potentija. i , which yields the effect of a finite nuaber of blades. 

A partial differential equation for is obtained and then converted 
into an ordinary differential equation by assuming a suitable form of 
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solutioa which is periodic in the circvnafereatial direction 9 , 

This ordinary linear differential equation is solved by well known 
KB approximation. 

The vorticity distribution representing the blade canber 
can be directly obtained fron the prescribed velocity distribution, 
so the first approximation of the camber line shape can be detenained. 
Using the $ solution, induced velocities due to the vortex distri- 
butionl^on the camber line can be calculated. At this point, it is 
possible to find the velocities induced on blade camber line, both 
along the camber line direction ? and normal to it in n direction. 
Clearly the sum of velocities in the 5 -direct ion due to mean flow 

and the induced flow should be equal to the mean of specified velocity 
distribution on upper and lower blade surfaces at each point. In 
general, this condition can not be satisfied unless a source distribution 
a ( 5 ) is also located on the camber line. This is the basis for 
finding 0 (c) . It can be shovm that the introduction of a(^') 
dianges the effective value cl Y(i;, which intura demands a rocal- 
culation of o(5). The method finds the converged values of o(5) 
and y(c) iteratively. Once these are known the blade thickness 
t(5) and. camber line shapes are readily found. Hence the blade 
prof Ho is completely defined. 

2.1 Meridional Botentiai How Solution for a Buct Tlow . 

loth the AnHjrsis and the Inverse problems n^d t^^^^ nerldicinal 
potential flow solution (without blades or blade forces) as an inpuv. 
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Ihis solution gi-ves a set of meridional atwamiines and the 
velocities along them. 

There are two methods in common tisage for finding the potential 
flow solution, (l) the relaxation or matrix throu^ flow method, and 
(2) the streamline currature method, 

!Che "relaxation or matrix throu^ flow method" consists in 
soi-ving the equations of notion in terms of the stream function. 

Mnit e difference techniques a3ce used fcsr the solution of the resxiLting 
partial differential equation. The flow region is covered by a fixed 
grid system formed by two families of strai^t lines or curves or a 
ccmbination. At each grid point the finite difference equation is 
satisfied, resulting in a set of simultaneous algebraic equations 
vfcich nay be solved using iterative techniques or direct matrix inversion. 

The alternative "Streamline curvatvire method" consists in 
solving the equations of motion in terms of several variables, the 
main one being the meridional potential flow velocity t In this 

method the entire passage is divided into a diosen number of stream- 
lines and normals to then. Then iteratively these streamlines and 
normals are modified by satisfying the velocity gradient eqi^it ion along 
the normals and continuity equation in each stream tube. 

Both methods solve the same flow equat ions hut differences ^ 
in te<iiniques introduce different operational constraints, and 
difficulties, Davis and Millar"*^ have given a systtoatic ccmparisoti 
of tl^se two different methods for the calculation of potential flow 
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from, the point of view of a tiser v/ho wishes to select one for 
deveiopuent or implementation. 

Significant effort towards optimising the damping factor and the 
cvnvatore calculations is required if the streamline curvature method 
is to be used operationally. itor the present work the competing 
matrix method was considered preferable. 

Ehe partial differential equation for incompressible, axisymmetric , 
potential duct flow in terms of stream function 'I' , is (Ref. 1g) 







( 2 . 1 ) 


jPor the solution of this equation boundary conditions are 
specified at the passage boundaries in terms of ij) . We desire to 
impose uniform flow far upstream and downstream, for which purpose 
the passage boundaries are suitably extended upstream and downstream, 
and the values of ij; associated with the uniform flow are taken. 

On the passage walls, of course, iij takes constant values. 


The version of the matrix method used in the present work involves 
the covering of the whole region with a fixed regular square grid 
system as shown in fig. ( la) , and writing a finite difference approxi- 
mation to the principal eqn. (2, l) at every interior grid point and 
grid points r®ar curved boundaries. This gives one algebraic equation 
for every interior or boundary grid pxaint in terms of xp values at 
that point and the nei^b curing pjoints or boundary values. G3ae 
resulting system of simultaneous equations is solwed by the G-uass— 
Seidel iterative scheme. 
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Ehe method offers fast coavergence, The scheme used yields 
second order accuracy and is inherently stable. This method is a 
powerful one and can probably be successfully adapted for ti^e blade- 
to-blade analysis problem. [[[The details of the above method are 
given in Appendix l[[} . 


2,2 Velocity Potential for the Perturbation Velocity 

as shoun in fi^.2(a; 

The flow throxx^ a blade rowVis assumed to be incompressible 
and inviscid. Since the flow is irrotatioual everyvdiere, the 
velocity potential. $ , for the absolute velocity satisfies the 
Laplace’s equation. In the orthogonal curvilinear coordinate system 
(o,B,0) shown in fig. 2(b), Laplace’s equation for ^ can be written 


3.3 i 






, + 1 ^ 

+ 0 + r ar „ 2 


ar 


az 


.■1 ii. + + 

^ aa ^ r a r .2 


1^ 
r a r 


ae 

^ , j 

az“ r^ ae^ 




= 0 


(2.2) 


here oi = constant are the streamlines of the duct flow and 0 =constant 
are equipotential lines of that flow. The defining equations of a , B 
are as follows, with S , S- the lengths along the respective lines 

da = r V ds (2.3a) 

: ^ r “ : 

' AS.': 'V'-- 

On rearranging teims of the above eqa, (2,2) and assuming that 
when the blade row is operating near design conditions, the velocity 
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u^, in a -direction can be neglected as compared to the velocity v„ 
in the p-direction, we get (Ref. 9) 


2 2 2 ^ " 

33 r v^^ 3 6 

Ihe potential 5 can be considered to be the sum of two parts, 
the velocity potential of the mean flow in the limiting case when the 
blade row consists of an infinitely large number of blades, and a 
perturbation velocity potential 0 , which yields the effect of finite 

number of blades, The former will independently satisfy eqn. (2.4) , 
hence the latter $ must satis:^ 



1 


2 2 
^ V 



(2.5) 


Considering the periodicity in 
expressed as 

0 = B( 3) e 


9 -direction, 0 can be 


( 2 , 6 ) 


where = mH, (m = 1 ,2,3,. ..) and ff is the number of blades. 

Substituting eqn. (2.6) into eqn. (2.5) , an ordinary differential 
equation for B is obtained 



2 2 


r V 


pjf 



(2.7) 


9 . 


Iccordirig to Eashiwabarar if tho nimber of blades U is 
relatively large, the ’»1KB method" may be used to integrate eqn. (2.7). 
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u^, in a-direction can be neglected as compared to the velocity 
in the P-direction, we get (Ref. 9) 


2 
3 $ 

2 

3S 


+ 


1 


sfo 

a 

p* 


0 


(2.4) 


Hie potential 5 can be considered to be the sum of two parts, 
the “velocity potential of the mean flow in the limiting case when the 


blade row consists of an infinitely large number of blades, and a 
perturbation velocity potential 0 , Vvhich yields the effect of finite 

number of blades. The former will independently satisfy eqn. (2.4), 
hence tie latter i must satis:^ 



2 2 
r T , 



(2.5) 


Considering the periodicity in e -direction, 0 can be 


eicpressed as 

0 = B(8) e (2.6) 

where X = mR, (m = 1,2,3,...) and is the number of blades, 
m 

Substituting eqn. (2.6) into eqn. (2.5), an ordimiy differential 
equation for B is obtained 



>1 

2 2 


r Y 





(2.7) 


n 

According to Ka^iwabara^ if the nxmber of blades IT 
relat ively large, the **VilKB method ’’ may be used to integrate eqn. (2.7). 
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The method gives an exact solution in the limiting case " • 

mien X is large, the accuracy of the solution will be sufficient 
m 

for practical purposes (Eef. 17). 

The WB solution to eqn. (2.7) can be expressed as: 


B = a^G^ exp (X^t) + exp (-X^t) 

( 2 . 8 ) 

where 


T =/p‘'/2 (jg 

(2.9) 

p = - 

r - 22 

( 2 . 10 ) 

"" V 


a , H axp Ci ^ 

( 2 . 11 ) 

A . = / dg 

( 2 . 12 ) 

dg 

(2.13) 

K 5 (4PP" - 5P* Vl6I^ 

(2.14) 


are constants. Subscripts 1 ,2 correspond to upper and 
lower signs respectively. Hereafter prime denotes differentiation 
with respect to 3 • Now the fundamental solution of eqn. (2.5) will 
become 


= ‘“mi 77. 

10 20 


(2.15) 


maere 


f (S) 


= C 


ml, 2 1,2 


exp ( 


+ X t) 

— m / 


(2.16) 
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The siibscript o denotes the value at the point (g^,0^) as 
shown in fig, 4f and » 0 ^^2 constants, Ercm. eqn, (2.15) 

the ^neral solution of 0 can be expressed as: 


m=1 m=1 10 20 


2,3 Induced ftelocities due to BLstribution of Singularities ELaoed 
on the camber Line of Blade Sections 


The introduction of a finite but sufficiently large nmber of 
blades is expected not to disturb the stream surfaces of the mean 
flow. In this situation the flows on the stream sheets can be 
considered two-dimensional, and the blade can be replaced by distributed 
sources,sinks and vortices placed on the camber line. The flows on 
the mean stream surfaces are however not independent of each other and 
we consider the relation between them in Section 2,6, 

In the above paragraph, mean flow surfaces refered to are the 
meridional stoream surfaces with an infinite number of blades present. 
These surfaces are not, however, known in advance. It is shown in 
section 2,6, that the difference between them and the meridional 
potential duct flow surfaces is small, imder certain conditions. These 
conditions are assumed to hold for the present analysis. Thus instea,d 
of analysing the three-dimensional flow passing through blades , we 
assume that the three-dimensional flow is composed of stacked two- 
dimensional flows, within the mean stream surfaces. 
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The potential of induced velocities due to sources, sinks and 
vortices placed on the comber line of blade sections as shown in 
fig, 4, is symmetric aboub 6 = 9^ with the periodicity of 2ir/5r , 
and tends to a constant value, say zero, as $->•_+ ® , Consequently 
from, eqn. (2.17) induced velocities can be obtained as given below, 

2, 3 (a) Induced Velocities in g and p directions 


Refering to fig, 3(a), induced velocities in C and n directions, 
vdaich are, respectively along and normal to the blade camber line, 
can be expressed as, (Ref, 9), |]]A condensed derivation of the 
following expressions is given in Appendix-II ]] 


C 

^ * 
P* 


a I ~ 
<' *^0 q V 

6i 


„ cos 
p*o 0 


(2.18) 
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^q *’ TTr_ V 


I lysine cos X„(0-0j + — ^ 4 X 


R f 
0 p*o m=1 0 o 


n 0 ^ rv „ f „ m 
0 


cos e sin 




-ia = 




/ T ^ 

'^o q V „ (x-x,) cos e, 


J H 

g irr^T 


■“‘o" 

00 - 

^ {x-xj 1 ^ 


o p*o 


"O' 3 f 

in=:l 0 o 


(2.19) 


'T~ Xm sin e sin X (0 -6„) } 
rv_^ f m o' 

p*- o 


C d3 

nv r , 0 


1 = r I . 

vp* ft . Y V „ / \ 

^1 p*o (x-x ) cos e. 


(2.20) 
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sin X (6-6^) 
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00 -pi P'1 
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Jy=~ir ^ X ^^^^"“^0^ 

m=1 10 20 o -^0 


( 2 . 21 ) 


^ |, Sin e oos X (6-e„)) 

p* 0 


Where x = N't , and f takes the form of f^ or according 

as 3 1 . °r 6 ^ 3^ , i.e. as 3 is upstream or downstream 

of a singularity. 

Now, if we define a new variable 0 by, 


X - X. 

X - x^ — - (l — cos 0 ) 


(2.22) 


then,eqns, (2.19) and (2.20) become Poisson integrals as given below, 


_|22. = I I j 

V 0 "1° 


Q IT 

=1 / I 


a r sin 

0 0 

9n 

d© 

0 


cos Sq 
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0^ - cos 

®c) 

ro ^o 

®o 



COS 

(^cos 

©Q— cos 



(2.23) 


(2.24) 


Hereafter si;hscfipt c denotes values at the c amber line. 
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•Ihe fianctions I , J , I , J in eqns, (2,18) to (2.2l) 

Q^C } C I C 

express the combined influence of all blade siogiiLarities , 

includicg those of the whole of the cascade. In the limiting case 

(b ->S*9 ->-0)» the influence of cascade vanishes, and the 
0 c^ o c ' 

terms J /(x-x ) and I /(x~x ) tend to infinity. The resulting 

0 C O I c o 

singular integrals can be evaluated for principal values only. To 
facilitate doing so nunerically a function T is introduced such 
that, 


where , 


lim 
f ^o 2c 
2o 9c 


_2£. 

x-x_ 


lim T 
I3c\ 

00 0c • 


lim 


8^ 9^ 
0 c 


x-x. 


lim 

Bo 


n ' 

Be 


V 6„ 9. 

^0 c 


T = 


ET cos e 


2^ Vp^o^^-’^o^ 


(2.25a) 


(2.25b) 


(2.26) 


It is easily seen that T expresses the influence of a 
singularity situated on the camber line, close to the singularity, 
i,e, in a region where the induced velocity due to the singularity 
is the deminant cbhtribution, 

Eqns. (2.18) to ( 2 . 21 ) give the mean -value of induced velocities 
on the caaber line (6 , 8 ). The induced velocities on both sides 
of camber line and close to it can be expressed as follows : 
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C “ Hh -10 ^ 

nqc ~ Q = + 1 2. + nqc 

p* p* p* 


(2.27) 




111 


0 * 


(2.28) 


Yflaere the + and - signs correspond to the suction and pressure 
sides of the camber line respectively, as shorn in fig. 3(a), 

2.4 Jeteimination of the Strength of Vortices and Source s/Sihk:s. 


Vortices ; 

A vortex sheet is equivalent to a discontinuity in the tangential 
velocity of a magnitude equal to the local strength of the voitex sheet, 
Ihis fact is used to find the vortex sheet strength y(C) along the 
camber line. Since the blade thickness is finite, thou^ small, some 
relationship is required to relate the surface speeds to the speed at 
the camber line of an infinitely thin blade. This is provided by 
the Eiegel’s relation (Eef. , 

The superscript denotes values relative to a blade and 
subscripts i ^ denote the upper and lower surface of the blade 
section as shovm in fig, 5 (a). voiticity distribution on the 

camber line is therefore, 

= t°5,o-0- * (1-7 (2.30) 
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The blade thickness t(5) is not known in advance, but values 
of t(C) depend upon cr(5), and are obtained from the previous 
iteration in an iterative solution process. 


Sotnrce/Sink Strength ; 

As already stated, these strengths are chosen to satisfy the 
condition that the mean surface velocities along the camber line 

have the same value when calculated in two ways, 

(i) as the mean of specified blade surface speed, corrected for 
the blade thickness by eqn, (2.29) i.e. 




(ii) as the combined velocity due to the mean and induced flows 


= v^ cos e + (w„-. /jor) sin e + C_ + C_ ^ 
5c * ^ ^yc 5qc 


(2.32) 


TAhere v.j^ and are the 0 and 6 components of the flow due 

to an infinite number of blades. 


Using eqn. (2.3b) and figs. 2(b) and 3(a) we can write as. 


dC = 


d0 


V ^ cos e 


(2.33) 


So the eqn, (2.19) can be rewritten as, 




I % 7 ^. 


(x~Xq) 


(2.34) 
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Using eqns, (2,25a) and (2,34), for convenience in evaluating 
0 , we write 


V 


£ J 


' “o “ ( “o < ^ '^0+ / o„ T dS 

p-)r 0 0^ 0 ''^0'^ DO c 


or 


^Q^O j U 


°Cqc,G 





T} dC^ 


(2.35) 


This can be seen to be, essentially, the induced velocity on 
the blade camber line, due to the members of the cascade excluding 
the blade concerned, or briefly the effect of the cascade. 
Similarly E is given by, 


— = / a ldC„ (2.36) 

V „ ■' o 0 

p* 0 

is seen to be, essentially, the induced velocity on the blade camber 
line due to singularities on itself. 


Now, using eqns, (2.31) to (2,36), E can be expressed as, 


_ 1 


rdt\ 2 / — 




) - y cose - (w„-cor) sin e 


“ ^5YC ” °?qc ,G 


(2.37) 


When the strength of sources and sinks 0(5) in the right 
hand side of eqn. (2.36) is assumed to be known, it becomes 



Fredholm's integral equation of first kind. Using eqns. (2.22) 
and (2,26) along with eqns. (2.29) and (2.33), the eqn. (2,36) can 
he rewritten as given in Eef. (9) , 


E 1 


V „ cos e 2 IT 

p* 



sin 0 d0 
0 0 

(cos 0^ - cos 0 y 


where 


V cose 

p* 


( 2 . 38 ) 


In order to get the closed blade profile, the following closure 
condition is to be satisfied, 

IT 

! I sin 0 a© = 0 (2.39) 

0 

The solution of a similar integral equation as eqn. (2,38), 
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which satisfies eqn. (2.39) has been obtained by Truckenbrodfc . 

The solution of above equiations is given on similar lines by 

9 

Kashiwabara as follows. 


I 


IT sin 0 


I 


B 


. 2 

sin 0, 


d0. 




(cos 0|^ - cos 0 ) 

here as E, given by eqn. (2.37), also contains the unknown 
0(5), so eqn, (2.40) must be solved by successive approximations. 


(2.40) 


2,5 Blade Profile Shape 

Blade profile shape can be defined by shape of the camber line, 
and the thidsness distribution along the camber line, Th^se can be 
obtained as given below. 
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2, 5 (a) Camber jjine Shape 

If the aagie of iaciden.ce at the, inlet to a blade row is zero, 
it can be expected that the row will exert negligible upstream 
infLueace on. the flow. Making this assumption, the velocity in 
O-direction, in the flow of an infinitely l:irge number of blades 
has the following relation with the circulation r(g) about the 
z-axis. 


r(5) = 2iTr(?)w^(5) (2.41) 

vdth r(5) defined by, 

r(5) = N / y(C) dC ( 2.42) 

0 

and an irrotational inlet flow is assumed. 


Assuming that the relative flow runs along the camber line and 
cannot penetrate throu^ it, i.e. the velocity in n -direction at 
the camber line is zero, we get 




-v^ sin e+ (w^-<or) cos e + 


(2.43) 


This eqi&tion will yield e , which can be used for obtaining 
the 6— coordinate of the camber line as follows: 

y p 


tan£= dSg/dS^ 


m 


e= / 


tan e 


r V. 


de + 9 




cx 


or 


(2.44) 
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Brom continuity equation between half the thidmess of the blade 
t and the distrib;xt:ion of sources and sinks, as shown in fig, 3(b), 


we get 


1 ^ t 

- / cr b d5 = / Cc bdn 


(2.45) 


where b is the distance between the mean stream surfaces. If t 
is assumed to be comparatively small, and using the velocity on the 
camber line, eqn. (2,45) can be approximated as follows, 


1 / abdS Oc bt 

2 0 “ 

where G is given by eqn. (2,32). 


(2.46) 


Since from fig. 2(b) , 


rv „ , hence 

p* ’ 


d a 

dS„ = b = 

a rv 


Becau^ da is constant between any two stream surfaces, eqn. (2,46) 


can be written as, 


1 / o : _ 

2 / r v„ =rv_„ ^ 


(2.47) 


This equation will give the thickness distribution along the 
camber line. In order to obtain closed blade profile the total 
strength of sources and sinks must be zero. Consequently the thickness 
at the trailing edge should be zefo. So from eqn. (2.47) we get. 
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£ 

/ d? = 0 (2.48) 

0 p* 

This is also a closure condition. The eqn. (2.39) is an 
alternative form of the above equation. 


2.6 An Estimation of v„ 

“?V 

It is shown in Appendix-III that if the blade force in 

9 

a -direction and (i^’w^) are both negligible, then v^= v . 

Jr 

The latter condition is seen to mean, on applying Euler's turbine 
equation along streamlines, that the work done on the fluid at 
different blade sections is so distributed that the total energy of 
the fluid along the 8=constant lines is constant. This condition 
can be achieved in ' free-vortex* axial flow designs. In general 
it is not achieved and v^ differs from v Eor small departxnes 

-V p* 

from the free-vortex condition, and with the ct -component of blade 

9 

force negligible, Kashi wabara gives the following relations: 


a 

ic= / 

a. 


_L 

2 

'p* 


(to - -“ ) (rw*) da + K 


h 


( 2 . 49 ) 


vhere k = 


(2.50) 


and is given by j 


w^ 


'h 


— f ! -j" ^ sT 

“h “h V 


da 


( 2 . 51 ) 



24 


where Q is the volvme flow aind has a relation, 

If =(0.^-0^) (2.52) 

!I!hus above e goat iore can be used for calculating using known 

values of and v^^ . derivation of above equations is 

given in Appendix-Ill 3 • 



chapter III 


PROCEPDEE FOR NUfffiRICAL SOLUTION OF THE INVERSE PROBLEM 

The proceduiB for the numerical solution of the inverse problem 
using a singularity method is given here. The analytical formulation and 
a brief outline of the procedure are already described in Chapter II. 

Before proceeding with the solution, many parameters are to be 
selected beforehand such as 5 geometry of the flmi passage of turbomachine 
i.e. of the hub and the casing? position of the inlet edge; design volume 
flow; number of blades and camber line lengths of blade sections on each 
stream surface. Velocity distributions are specified on both blade surfa< 
as a function of camber line length. From the generalisation view point, 
all values are used in the non-dimensional form with respect to the 
inlet radius at the casing and blade speed at this position. For the 
numerical example presented here values of above parameters are given 
in Appendix-V . 

In Fig. 5 the flow passage of a mixed flow turbomachine, five 
stream surfaces of axisymmetric potential duct flow and the blade inlet 
edge are shown, which are cited from Ref . (9) • la the numerical example 
solved here , blade geonetries are calculated with the prescribed velocity 
distributions show* in fig. 6 , also cited from Ref .( 9 ) , for all the five 
stream surfaces consiife red. 

The stepwise procedure for the numerical solution of the inverse 
problem, based on the theory given in Chapter II is as follows 
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(1) The meridional potential duct flovv solution without blades or blade 
forces for the given turbomachine is found using the method chosen 
in Section 2.1. This will give a set of streamlines and velocities 
along them, which are used as an input to the next part* The details 
of the procedure for finding the meridional potential flow solution 
using a matrix method are given in Section 2.1 and Appendix-I* 

(2) Lengths of streamlines are calculated from the coordinates of 
streamlines, starting from the given position of the inlet edge. 

(5) As is well known, a good specification of the hlade geometry should 
contain relatively more data for the highly curved nose and tail 
regions than for the rest of the blade* Further, in the interests 
of limiting computer storage and running time demands, the blade 
geometry is sought to be specified at’ the^ fewest possible points — 
nineteen in the present case. Glauert's transformation = 

(1 - cos <f>)/2 yields a suitable distribution of points in physical 
space, when eq,ual intervals are chosen for . The given table of 
values for the velocity distribution is transformed to 

table of values at the chosen <f> points by interpolation. 

(4) Prom this specified velocity distribution along blade surfaces the 
first approximation of vortex distribution along camber lines is 
calculated from eqn. (2 .50) , assuming zero blade thickness . As 
the circulation of the flow at inlet is zero, from e<in.s. 

and (2.42) % can be calculated. 

(5) How assuming that v^ snd v^^ are the same for the first approximati 
eqn. (245) is used for cal onlating the angle e of the camber line 
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with S ^-direction, neglecting induced velocities . due to vortices, 

u* 

sources and sinks. The camber line slope tan e , will be used to get 
the coordinate Sg of the camber line. The eq,n. (2.44) will give the 
6 -coordinate of the oarber line at corresponding points. These S 

p 

and ©-coordinates are the first approximation to the shape of camber 
lines. These camber lines are taken as relative streamlines on 
respective stream surfaces . 

(6) Variables appearing in the WKB solution t, P, K, f^ and f^ are 

calculated using eqns. (2.9) to (2.I4) and (2.I6) . Derivatives of 

these variables P' , P' , K' , f and f' are calculated using suitable 

1 2 

numerical techniques. 

(7) Using the vortex distribution obtained above and variables 

along with their derivatives computed in (6) , the induced velocities 
due to vortices placed on the camber line are calculated from eqns. 
(2.21) and (2.24) • The infinite series appearing in these equations 
should be truncated suitably. 

(s) The prescribed velocity distribution and eqn. (2.37) are u.sed 

together to solve the integral rc u. (2.4O) . In the integral eqn-(2 
the unknown 0(5) appears on both sides. This integral equation mus' 
therefore be solved by successive approximations. The solutions bn 
different stream surfaces yield the source-sink distribution cf(C) 
along the camber lines on respective stream surfaces 
(9) Induced velocities due to source-sink distribution <7(5) are > 
calculated using eqns. (2. 18) and (2.23) • Here also the infinite 
series appearing in these equations should be truncated, suitably . 


( 10 ) 


2B 


Now the shape of the camber line is corrected according to the 
condition that the relative flow runs along the camber line. 

Eq.ns. (2.41) to (2.44) along with eqn. (2.30) are used for this 
purpose. This step is quite similar to step (5) . 

(11) Prom the source-sink distribution obtained in step (s) , the blade 
profile thickness is calculated using eqhs. (2.52) and (2.47)* 

(12) N07ir v^ is evaluated from knwm values of v^^ and on each stream 
surface, as given in section 2.6, using eqns. (2.49) to (2.52) . 

(13) The vortex distribution is calculated again from eqn. (2 .50) taking 
into account the blade thickness obtained in (II) . 

The stepwise procedure given ^ ove from step (6) to step (13) 
is iterated until the shape of the camber line and blade thickness 
converge with desired accuracy. 

It is to be noted that step (12) given above is included in the 
iteration loop. The reason is, that w* will change due to successive 
approximation of blade thickness. Her.ce the iterative calculation of 
is also needed, because eqns. (2.49) (2.5’') contain as one of the 

variable*. 

The procedure described above will remain the same whether the 
hydraulic turbcmachine is a turbine or a pump, except' that the angular 
velocity term in tbe analytical fcrmulation will be -o) for a pump. 


CHiPTER I\r 


mERICiL lECMIQUES TJSED AND 13®]]?. SUITiBILITY 

The adoption of numerical techniq.ues for the solution of engineerin 
prohlems has got an important place, hut considerable attention must be 
paid while using these techniques. 

The use of numerical analysin is both a science and an art. It 
is a key to the specialist in the field but is often misunderstood by 
nonspecialists. Like most arts it is learned by practice. Principles 
there are, but even these remain unreal until actually applied. 

A person wishing t* solve an equation or some other problem 
numerically, usually has a wide choice of methods. Most problems of 
engineering interests are not so generously endowed, and for them the 
choice is restricted to a few methods. A reasonable balance must be kept 
between memory or storage requirements, accuracy, and speed of computati( 

In the present numerical solution, the above considerations were 
kept in mind while selecting any particular numerical technique, by the 
limited knowledge gained by the author. The numerical techniques used 
along with their suitability for the present numerical solution are 
discussed below. 

4 .1 Solution of. Simultaneous Linear Algebraic Equations s 

In the meridional potential duct flov/ solution, a system of 
simultaneous linear algebraic equations results from the finite differei 
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approximation of the partial differential eqn. (2.I) at each grid point. 
The number of grid points is very large (about 2000 for the present case) , 
for getting an accurate solution of eqn. (2.I) in the given flow region. 

It is almost impossible to store a coefficient matrix of this size 
(2000 X 2000) in any computer. 


One important property of such a system of algebraic equations 
is that it has, at the most, only five non-zero terms in each equation, 
and the coefficient matrix contains a high proportion of zero elements. 
Consequently attempts to work entirely in terms of matrices and, in 
particular, to store all their elements is either impossible or extremely 
wasteful of computer storage. Bather, here the Gauss-Seidel iterative 
method can be iised by the repeated application of a single algebraic 
equation throughout the whole system of grid points . Using double 
space subscripts are shown in-fig. 1 (a) , this amounts to computing a 
new approximation of . at every interior grid point from a simple 
formula. For the present case, e.g.. 


+ii + ''P . ('P . . ^ . ..) Ax 

^ _■ i-1, .l ^1+1 ,.a±l _ ^ 

^i,3(new) 4 

... (4 .l) 


where Ax is the mesh size. For grid points near the boundaries 
similar equations are used, which are given in Appendix-I, along with 
other details. 

4,l(a) Acceleration of Convergence s 

For accelerating the Convergence of the iterative process used 
for the solution of the system of algebraic equations arising from the 
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elliptic problem, the method of successive overrelaxation is used# It 
is again an iterative technique, in which an improved estimate of 4). . 

is computed by applying the following formula at every grid point* 

h, 3 (improved) 

where, is an arbitrary parameter. For R^ = 1 , the process is identical 
to the Gauss-Seidel method. However, for a choice of R^ in -the range 
of 1 < R^ <2, the convergence is more rapid for elliptic problems. 

Finding an optimum value of the parameter, , is somewhat difficult, 

20 

and is discussed in detail by Forsythe and Wasow . For the present 
case it was found by trial, and its optimum value is approximately 1.8. 

It was observed that number of iterations req,uired for a chosen degree 
of convergence falls very rapidly, when the parameter is in the immediate 
vicinity of R^ , and it is generally better to over-estimate Qp^ 
than to under-estimate it. 

4 .2 Interpolation s 

In the meridional potential duct flow solution, a second order 
interpolation is employed near the boundary grid points for evaluating 
the functional values near the boundaries. As there we want functional 
values in terms of the first and second partial derivatives at the boundary 
grid point, so a Taylor series around the boundary grid point is a suitable 
interpolating formula. (Ref . ii) • 

Jhrerywhere else, for interpolation in terms of only functional values, 
a Lagrangian q,uadrati6 polynomial is used, because the points are uneq_ually 
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spaced* Acton suggests that the Aitken' s interpolation is (luite 
suitable for unequally spaced points. It is no doubt superior in the 
sense that in Lagrangian interpolation there is no indication of the 
accuracy, whereas in Aitken' s interpolation there is. However the 
accuracy of Lagrangian interpolation differs usually in 5th or 6th 
significant digit only as compared to iitken's interpolation. Therefore 
the Lagrangian quadratic polynomial is used here, because of its simplicity 

4 .3 Differentiation ; 

The calculation of derivatives numerically is a hazardous 
operation, especially when dealing v/ith low accuracy data. And even for 
high accuracy data, derivatives higher than first are likely to have 
considerable errors. The inherent difficulty in numerical differentiation 
is that it tends to magnify small discrepancies or errors in the 
approximating function used for representing the data. That is why 
numerical differentiation is often called a noise-magnification process 
(in contrast to integration which is called a smoothing process) . Because 
of the tendency to magnify discrepancies or errors,- numerical differentia- 
tion is usually avoided wherever possible. 

If it must be used considerable care is required if serious errors 
are to be avoided. 

A number of methods are availaole if the data is equally spaced. 

And for many methods it is quite possible to estimate an error bound y 
for the derivative computed* But for unequally spaced data it is quite 
difficult, if net impoasibie ,to get an idea; of the error -present^ rte 


derivatives 
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In the present numerical solution) almost everywhere the data is 
unequally spaced, and the derivatives even of the order four are needed. 
Keeping this in view, the selection of the differentiating method was 
very carefully studied. It was foiind that no available numerical method 
can control sizable errors in higher order derivatives. So a variety 
of formulae were tested on a sine cum^e to select one with less error in. 
the higher order derivatives, along with reasonable storage requirements 
and speed of computation. It was necessary to impose the last two constr- 
aints as the derivatives are to be calculated frequently, and the solution 
as a whole needs a lot of storage. 

22 , 

Wilkinson made an attempt to compare many formulae for 
calculating derivatives numesrically, particularly for second order 

derivatives, using equally spaced data, ’.^rhere he used a cosine curve for 

c 

the above purpose. For unequally spaced data Wilkinson suggested a method 
of "parametric differentiation". In this method, both the dependent 
variable y and the independent variable x are assumed to be functions of 
the point index number i, and then formulae for equally spaced points are 

used to differentiate both the variables w.r.t i. Formulae suggested by 

S ' 

Wilkinson are as given below : 




dx d^y dy d^x 
(dx/di)^ 


(4.5) 

(44) 


Even this 'parametric differentiation method* , when used for third 
or higher order derivatives, was found to result in sizable errors in the 
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derivatives of the sine curve used for testing purposes. It is expected 
that if the variation in the spacing of data ,is small and smooth, the 
results will be somewhat better, as the errors appearing may be due to 
the lack of smoothness of variables with respect to i. 

It v;as fomd that fitting a second order polynomial through three 
points to find the derivative at the mid-point is somewhat superior 
(i.e. the error incurred is comparatively smaller) . As in the case 
of equally spaced data when second order polynomials are used, the error 
in derivatives at the endpoints is almost twice that at the midpoints. 

It was found that a simple two point linear differentiation formula for 
the end points gave a reduced error as compared to second order polynomial. 
As this combination gives better results for the derivatives of an 
analytically knovm function (sine curve) , and since the st:)rage requirements 
and speed of computation are reasonable, it was used throughout for the 
calculation of derivatives in the pr-sent case , Derivatives higher than 
first are found by successive differentiation. 

The formulae used for computing derivatives are given beloi,?. These 
are very simple and can be derived very easily. 

At the end points, 

(^1 “ (j-a - 

For the las t jjpint formula w ill be s imilar to (4 .5) . 

For rest of the points, 


(4.5) 


(4.6) 
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where , 


^M-1 _ %+l 

hi hi2 hi h^ ^ ^^2 


( 


(4.7) 



=ij+i^ 


*^1 ’'12 



>^M+1 




^2 N 2 


n 


... (4.8) 


and h, , hj - 


V N2 ■ ^ ■"‘'2- 


4 .4 Smoothing : 

The meridional potential duct flow solution gives the streamline 
coordinates (z,r) and meridional velocities v^ along them. As in the 
present case, first and higher order derivatives of K and P (defined by 
eqns. (2*10) and (2.I4) , which are basically functions of r,z and v J 

ir 

will be used in subsequent calculations, it is better to use least square 
polynomials to smooth the data (r,z,v , before calculating K and P and 
differentiating them. 

s. ■■ 

A five point smoothing formula adopted by ITilkinson for smoothing 
the value at the midpoint is, 

yjg == _ 0.00174 yM„2 + 0*^875 y^^i + 0.62847 + 0.1875 yM+l"® ^M+2 

... (4.9) 

This formula is used for smoothing r,z and v^^ for each set 
except two points on each side of the range, which remain unchanged. : 

4 *5 Integration t 

is unequally spaced almost everywhere , so a composite 
integration fornnila, using a second degree interpolating polynomial passing' 
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through three consecutive points is employed for integration on sub- 
intervals. The repeated application of a low-order formula is usually 
preferred to the single application of a high-order formula, partly because 
of the simplicity of the lovir-order formulas and partly because of 
computational difficulties. 

The formula used degenerates to Simpson's rule, for equally spaced 
data. Its derivation is quite simple and straight-forward. 

4 *6 Fourier Approximation, Discrete Domain s 

Poisson integrals appearing in eqns. (2 .25)? (2.24) and (2 40) 
are singular integrals, and the function inside the integrals is highly 
non-linear. Because of presence of singularity these integrals must be 
handled carefully. 

To evaluate these integrals numerically as given in Chapter V, 
the discrete Fourier approximation ■ is need to represent the integrand. 

For the present case, only the Foijrier sine-series approximation for the 
half range is sufficient , For the discrete domain 0 to n, the formulae 
given below are based on Ref. (23,24) . The function g (0) can be 
expressed as, 


W 

n=1 


(4.10) 


M = 1,2, .. .,N+1 


h-l 


\ 'h Sinn 

1=1 


(4.11) 


n ^ 1; 2. j • ft • « 
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where (N+1) is the number of points at which function is known, HI'’ is 
the number of terms of the Fourier series, less than H. Using (H-1) 
terms of the above Fourier series (4*10) , the function can be represented 
exactly at all (H+i) points. But if smoothing is desired, one can use 
less number of terms , 

4 .7 Shape of the Camber Line ; 

7/hile determining the shape of the camber line, the angle € is 
to be calculated using eqn. (2.43) > v/hich is a transcendental eq.uation 
for e . This equation could be transformed rather easily to the form s 

p 

sin e + B sine + C = 0 (4»‘12) 

where B,C are known functions of (‘B, 6). Thus e can easily be found. 


However an alternative approach making use of the perturbation 
principle, as given below was also used. 

As the induced velocities C and C are small compared to 

pYc nqc 

other terms of the eqn. (2.43). So an initial approximation e^ can be 


obtained as, 


w -mr 

s, = tan-^ (-— ) 


(4 .13) 
(4.14) 


Defining ^ ” V ^ 
and using the usual perturbation procedure to find 5 from eqn. (2.43) » 
one gets 


C + U + (w^ - (ijr) :cbs - V* sin 
DYc nqo * ^ o ^ Q 


v^ cos e^+(w^ -m r) sin 


( 4 . 15 ) 
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The perturhation process can he applied repeatedly to get a 
converged value for o • However, it v/as found that only a single shot 
is sufficient to give a good approxiina.tion for e « 

4-8 Length of Streamlines in feridional Plane s 


Length of each streamline in meridional plane is to he calculated, 
as it is used for the solution. It is calculated from the leading edge, 
using the formula, 



where (r,z) are coordinates of a streanaine, and dr/dz is found by fitting 
a local quadratic through three points ( D » l+'l > j+2) , which gives 
r = Az + Bz + C, and dr/dz = 2Az+B. The eqn. (4,16) gives 

^6 = ^8 + dH (4.17) 

Vi H 

J 

where H 2Az + B, hence 


^0+1 


^6 ^ 4A 
J 


QH '4 + 


+ Aog^ (H + 


'4 + H^)“I 


D+1 

H. 

.1 


(4.18) 


If the value of A is ver5'' small i.e. the streamline is locally a 
straight line, the following formula will be used instead cf eqn. (4 .18) , 
which is obtained from eqn, (4 .16) as, 



CHAPTER Y 


COMPUTATIOITAI, DIFFICULTIES AIID THEIR EJMEH- 


As indicated in Chapter lY, the use of numerical analysis is 
both science and an, art. The science and art juxtaposition is due 
to an uncertainty principle which often occurs in solving problems, 
namely, that to determine the best way to solve a problem may require 
the solution of the problem itself. 

The art aspect of solving problems numerically arises mainly 
in two places, (i) in choosing the proper method, and (ii) in circum- 
venting the main road-blocks that alw^s seen to appear* This 
implies the need for anyone who wishes to use numerical techniques to 
deveiope experience and with it, hopefuCLly, intuition. 

Some of the computational difficulties encountered during the 
present numerical solution along with the techniques used to circumvent 
them are described below. 

5.1 Evaluation of Poisson Integrals 

Ihe Poisson integrals appearing in eqns, (2.23), (2.24) and 
(2,40) are singular integrals. Also the function inside the integr^s 
is hi^ly non— linear. For evaluating these integrals, proper care 

is needed because of the presence of the singularity. The procedure 
adopted for this purpose is as follows : 

The Poisson integral can h® written as, 
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/ g( 


0 



sin 0 d 0 
o 

(cos 0^ - cos 0 ) 


(5.1) 


For the present case, the function g ( 0^) is known at 

nineteen points, uneqtaally spaced. Interpolation is used to get 

values of g at intervals of n/l8 , that is at nineteen equally 
spaced points. At both ends of the range of integration the value 
of s( Qq) is zero {which is obvious from the form of g (0^) in 
eqns. (2.25), (2.24) and (2,40)} . It can therefore be represented 
by a half~range Fourier sine-series, in the discrete domain. The 
Fourier approximation of the function g( ), over the discrete 
domain, 0 to u , is given by eqns. (4.10) and (4 .II). Using these 
equations, the integral (S.l) will have the form, 

ir NF sin 0 d© 

/ { y B sin n 0 } — (5.2) 

o n=1 (cos 0 - cos 0 ) 

o 

In this integral the 0^' constants, so it can be 

written as; 


ITF 

y { B / 

L n ■' 
a=1 


TT sin n 0^. sin 0^ 


0 (cos 0Q ~ cos B ) 


d6„ } 


(5.3) 


with 


sin n 0 sin e = -1 { cos j^Cn-l) ® “1-cos l"(n+l) 0 yiy^t^ 


The integral becomes a Sum of integrals of the foim, 


TT cos n 0. 


o (cos 0^ cos 0 ) 


. d0. 


(5.0) 
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Now, this type of integral may be eyaluated using the 
following resiO-t, derived by Glauert and given in Ref. (is), 


^n= / 


IT 


COS n R 

o , sin n 0 

— — de =1T— : — r 

0 (cos 0 - cos 0 ) ° ® 


( 5 . 6 ) 


Ihe recurrence relation of such integrals, given in 
Ref. ( 16) is, 


I _ = 2 I cos 0 
n+1 n 


- I 


rt-l 


( 5 . 7 ) 


where I^= 0 and = ir , this is obvious from eqn. (5.6). Ihe 
constant coefficients appearing in (5.5) can be evaluated 
using eqn. (4.I1). 


5.2 Convergence of Solution for the Integral Equation 

The integral eqn. (2.40) has the unknown variable a , both on 
the left hand side as well as inside the integral in the form of 
variable It must, theref ore ,be solved by successive approjdmations, 

which is an iterative process, A single iteration for the solution 
of this integral eqi^at ion at one stream surface, takes about 4.5 seconds 
of computational time for the ccmpvtter used, which is quite significant 
if a large number of iterations is needed to get a converged solution. 

So obviously, one is interested in accelerating the convergence of 
the solution^ The tricks employed for doing so in the present work 
are described 

The method of successive over relaxation for accelerating the 
convergence was found quite suitable. The formula for over-relaxaticm 
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given by eqn. (4.2) was used but the parameter i.e. relaxation 
factor, was not kept constant, as for the meridional potential duct 
flow solution* Por getting optimiJini convergence rate was varied 

from 1,55 to 1.75 (found by trial), depending upon the percentage 
change in the value of o , fccm iteration to iteration at any point; 
i.e, if the change is large, E^ will be low or vice-versa. 

With successive over-relaxation the convergence is quite fast 
for some but not all stream surfaces, if we start with initial approxi- 
mation as zero source distribution. In particular much trouble was 
encountered with the hub streamline, which also happens to be the 
longest. It never became clear, inspite of much experimentation, vtoy 
this streamline was so troubiescme. There was, however, some evidence 
that convergence was more rapid on the shorter streamlines i.e, those 
for which the mesh points v/ere relatively closely spaced. The 
procedure finally adopted was based on the expectation that the pattern 
of source/ sink distribution would be the same for all stream surfaces, 
and that an initial distribution displaying such a pattern would be 
beneficial. So the source distribution of the fifth stream surface 
was obtained first and iised as an initial approximation for the fourth 
stream surface and so on. It was observed that the convergence for 
the other stream surfaces was improved considerably using this proceduie 
together with the successive over— relaxation discussed above. 

It was found that the convergence can be improved further, 
if we use successive over— relaxation after three unaccelerated itera- 
tions to allow the solution to "stabilize” itself. 
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5.3 jBloating Point Overflow/Underflow 


Ydiiie evaluating jf2 cally, using eqn. (2. 16), the 

exponents of the exponential term of the order of + 200 and above 
appear, For the computer used for computation, the limit on the 
values of this is about + 88. 


On carefully looking at the equations, it becomes clear that 
whenever f^ or fg or their derived variables appear, they are in 
quotient form and the expected valves of such quotients are reasonable. 
How fcom. the computational view point eqn. (2. 16) has the following 
form. 


^1,2 “ ^ v5.8) 

hence f’ = (aBI, , + A*) exp (b J (5.9) 

Using eqns. (5.8) and (5.9), the different types of quotients 
to be used for the computation can be expressed as. 


^ exp 


f 


B’ + a'/a 




A + A!. ^ 

r ^ exp (B^^) 


fT L + AV 


(T> 


(5.10) 

(5.11) 

(5.12) 

(5.13) 
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^ A' I1«P 

=^0 - 


(5.14) 

(5.15) 


In eqns. (5.10) to (5.14) subscripts 1 or 2 of f, f ’ , 
f » , B, B’ , and B^ are omitted, as abo-ve eqi;&tions are true 
for both subscripts. 


While -using the abo-ve equations only the -values of B^jBg, 

B* - B’ . A and A’ are to be stored, and eqns. (5.10) to (5.15) are 

Y 2 * 

-used whene-ver appropriate quotients occur. 


Using this procedure it was found that the exponent tena in 
the quotients is always -"ve and is reasonably small, which makes the 
computation simple and reliable. 



CHAPTirp. VI 


EESDITS AITP DISCUSSION 

The americal solution for the inverse problem given here is 
for a mixed flow turbomachine shown in fig. 5. The computer program 
prepared is explained in' detail in Appendix - IV. Although written 
for a particular geometry, it can be readily modified to deal with 
any type of hydraulic turbomachir®. 

The turbomachine , velocity distributions on all five stream, 
suirfaces and other data are cited from Ref. (9). In fig. 7, the camber 
line shapes obtained are ccmpaied with those given in Ref. (9). 

Slight disagreement is present. It should be mentioned that the camber 
line shapes obtained had "almost converged" . Dm® to time limitation, 
the program was not run to achieve convergence within the specified to ler 
The blade thicknesses obtained by the present numerical solution 
are compared with those given in Ref. (9), as shown in figs. 8 and 9. 
There are considerable diffe3?ences in magnitude, but the trend of 
thickness distribution is similar. Much effcmt was spent in trying to 
locate the cause of this deviation, as described below, but without 
success. The iteration to iteration variation of blade thicknesses is 
given in figs. 8 and 9. It can be seen tMt thicknesses displayed 
a converge nt trend. linB limitations rather than any program defect 

prevented the solution being carried to a hi^r degree of convergence. 

The deviation in the results as compared to those given in 
Ref.(9) , may be due to some of given below : 
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(1) Jbcp©i'XD}Pn.'fca’tion. showed ■tha'fc "tli© solution is yory sons it iv© to 
yaluos of input data, particularly that of voiocity distributions* 

As an example of this, a slightly varied velocity distribution was 
input for the first stream surface. The results of camber line 
shape and thickness distribution for the first stream surface are 
given in fig, 10 for the two slightly different ^^locity distribu^ 
tions. The changes are notable. Clearly a precise agreement with 
the results of Ref, (9) could not be expected since input data 
were read off rather small graphs given in Ref, (9). 

(2) Further, it was found that even 5 to 10^ error in meridional 
potential duct flow velocity is sufficient to make the 
solution non-convergent, v is one of the important basic inputs, 

Jr 

and is not given in Ref, (9), The meridional potential duct flow 
solution used as input was found as described in' Appendix-I. 

There is reason to believe that on some points near the botmdaries^ 
values may change by 5 to 1^0 depending upon the detailed 
numerical method used for the solution. This could be another 
factor c arising disagreement. 

(3) Btoi* the results given in Ref. (9) , Kashiwabara^ used upto 30 terms 
in the series for the calculation of the perturbation velocities, 
whereas for the present numerical solution only 11 terms are iased, 
because of computer storage limitations. 

(4) Kashiwabara" has used 4 terms in the series for 0^ ^ ^ eqn,(2,ii). 
The third and fourth terms involve calculation of third and fourth 
derivatives of P, which can have sizable errors, and can make the 
soluction meanin^ess. , The numerical- methods examined for evaluating 
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these higher derivatives gave unacceptable results. So in the present 
work only 2 terms in the i®CB solution series were used. This difference 
in the number of terms used must clearly have affected the solution 
obtained, 

(5) For the estimation of v^. from v ^ using eqn,(2.49), the integration is 
to be carried out along 6 = constant line in the meridional plane. In 
the present numerical example it is assumed that the points having the 
same index number at every stream surface are on the constant line. 
This assumption is not perfectly valid, which is clear from the fig, 5, 
whe 3 re the leading edge is seen to be g = constant line, but not the 
trailing edge. Because numerically it is rather cumbersome to locate 
the 6 = constant lines and then evaluate the integral, so for simplicity, 
the assumption stated above was used with the expectation that .the change 
in v^ from v^^ will be small. This assumption can affect the numerical 
solution to some extent. 

Conclusions and Scope for Further Work ; 

The present method for 1he numerical solution of the inverse 
problem, was selected with the expectation that the method will be 
very useful. But the method as programmed in the present work, has 
tuirned out to posses some shortcomings. It is not however clear that it 
is inferior to other existing methods for the solution of the inverse 
problem. 

The effect of incorporating tbe corrections indicated in the 

foi^ going paragraphs needs to be studied v®ry closely. There is also the 
need to continue seeking better methods for the solution of the inverse 
problem, efforts for the same daould be continued. 
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IPOTDIX - I 


matrix method for 'THE IvERIDIONAL POTEHTIAL DUCT FLOW SOLUTION 


Tlie matrix method involves covering of the region, of interest 
with a fixed regular grid and writing a finite difference approximation 
to the principal eq_n. (2.l) at every grid point as given below. 


1 - Finite~dif fere nee Approximation : 


25 

Marsh used a curvilinear grid, because for the flow through 
a turbomachine of arbitrary shape, a quite suitable form of grid is a 
distorted or non-rectangular mesh, which is straight in one direction. 

He pointed out that v;hile using cuiwilinear grid there is no difficulty 

in forming the finite difference approximation for derivatives in the 

\ 

direction of equally spaced points, but there is no general expression 
for the derivatives in the other direction and these are to be found 
using Taylor series approximation to f^t a required degree of accuracy. 


In the present numerical solution square grui is used since this leads 
to simple expressions for the interdependence of the functional values at 
neighbouring grid points. 

At an interior grid point, using the notations shown in f ig. 1 (a) , 
the finite difference approxiBBtion of principal eqn. (2.I) is given by. 




^i-1,i ■ ^i-tl .d ^1, 3-1 ^1, 3+1 _ . 

4 :8r 


( 1 ) 


where ax is the mesh size. 
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Not/ for boundary grid points, consider a square grid shown 
in fig# 1 (c) , and let the boundary conditions be of the Dirichlet type. 
That is 4' is kno',m along the curved boundary, of which we shall use points 
A and B, For simplicity, abandoning subscript notation temporarily, 
the relevant grid points are labelled as 1,2 and 5, shown in fig* 1 (c) . 
For finding finite-difference approximation of eqn. (2.1) an appropriate 
Taylor' s series expansion is used and unv/anted derivatives are eliminated. 



B 


" * 3z 






( 2 ) 

( 5 ) 

( 4 ) 

( 5 ) 


2 ^ 2^ 

These equations enable us to find ^ and . The 

principal eqn. (2#l) then giwes : 


3r 
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= C ( D 'p 2 ® ^ \ + <> ’f'p) 


where , G = 


F = 
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1 _b__^ IT - 

» : ? = ^ 2r:t1 +by ’ 


a-tb 

ab 


(1-b) Ax 
2b r; 


1 

1+a 


1 - 


a 


(t+^T ’ ^ " b(1+^ 


Ax 

2r 


(6) 
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The method adopted for boundary grid points is a second order 
interpolation using Taylor’s series. It was found that the above method 
is quite accurate for the cases analysed. 

2 . Grid Size and Shape ; 

While using curvilinear grid, Marsh^^ pointed out that it is 
difficult to give a general formula for the dependence of truncation 
error on the local shape of the grid. However this error will increase 
as the grid is distorted from a squai^e shape. Similar is the case here, 
since the grid is not of square shape near the boundaries. Hence for 
such cases it is advisable to reduce the mesh size as the grid distortion 
increases. This will increase the amount of computer storage and 
computer time needed if a desired degree of accuracy is to be achieved. 

3 . Types of Boundary Points ; 

As shown in fig. 1 (b) , each point in row j can be placed into 
one of tte following four categories s 

Type I ! Point (i ^ l) , if this is not of type III or IT. 

Type II s Points (i, j), if any, for which > i > 

Type III ; Point (i i) , only if it is not of type IT and if 

^max,j+1 ^max,3 * 

Type IT s Occurs if the curved boundary produces an intercept within 

+ 6 of a grid point, where 5 < < 1. 

Hote that for any rcw j,, points of types II, HI and IT are 
mutually exoulsive . 



appendix - II 


ETIUCED YElOCITIES due to SUGDLAEITlES PLACED ® CAMBER DINE 


1 • Induced Velocities due to Source s/Sinks; 


Velocity potential of the induced Telocity due to a row of 
sources or sinks placed on the camber lines of N blade sections, as 
shown in fig, 4 , is an even function about 9=0^, with the periodicity 
of Eir/N, and tends to zero as 3 Consequently from eqn. ( 2 . 17 ) 

it can be exptressed as follows : 


ni=1 ^10 


« = I ■:?- cos X„ ( 6 - e J if 6 < B 




m 


00 

- y c 

" ^ m2 f^ m ' 0 

in=1 20 


(1) 


cos X ( 0-6 ) if B 1 B 


Where subscript q is used for ralues concerning soiirces/siriks, and 
subscripts 1,2 denote respectively the values upstream and downstream of ■ 
a singularity. As B= B^ is the matching plane between the upstreaa and 
downstream flow fields, the following continuity relations must be satisfied. 

At B= Bq (6/ 

; : :^q10 = ^20 ; : 

r (t - T ^) = o s Cb - ej 

©*^^20 qio ^ ^ 0 

where q is the strength of source or sink, and 6 is the Dirac s delta 

Equation ( 3 ) embodies the velocity discontinuity at the source. 


f met ion. 
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CDhe itourier expansion of a row of 6 -functions can be obtained as, 


— T V / \ 

Ce-SoII^ I (0-8 ) 

m=1 


( 4 ) 


Differentiating w.r.t. 0 will give rw^. How using the condition given 
by eqn. ( 2 ), we ^t 

0 In - 0 — 

*10 " ^20 

f f 
1 2 

as at 6 = 6 , — — = -ts- = 1 , we get 
^10 ^20 


"" °m2 


( 5 ) 


Similarly differentiating $ w.r.t. 6 will give y and then using the 

Q. Q. 

condition given by eqn. ( 3 ) along with eqns. ( 4 ) and ( 5 ) , Kashiwabara^ gets 


where P = 


Sal 

= "m2 = 7 

i 

fU 


20 

.lo 

^23 

^10 


W 1 


o ’'p»o ’ 


( 6 ) 


An. aticmaly in Kashiwabara^ s foimuLation of th6 problem is to be noted here* 
Y as obtained from (i) is a sine series in 0 / with a zero mean value 
when averaged for constant 6 * A row of sources on the other hand 
neces sarily gi^ves rise to a flow with a non-zero mean value j as is 
evident from the first term in(4)y aad thus can not be represented by 
a potential of the foim ( l)# Kashiwabara^ s formulation can probably be 
corrected in this respect. However, the Ixmited initial objective of 
the present work was to obtain agreement with Kashiwabara^s results 
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and this anomaly was ignored. Equatioos (5) to (s) are exactly those 
given in Eef. (9) ,where their derivation is not given. 


Using constants 0^2 gi'vea by eqn. (6), we get 
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V ^ 0 Vo m=1 ^0 V “ ° ° 


f. 


5 o Vo 111=1 0 -^10 
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( 8 ) 


where § is the length along the camber line from the leading edge as shown 
in fig. 3(a). 

2. Induced Velocity due to Vortices : 

The potential of the indiced velocity due to a row of vortices 
placed on the camber lines at shown in fig. 4, is an odd function about 
6 = 0^ and tends to zero as 6-v+<», so from eqn. (2.17) it can be expressed 
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as follows 
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\yhere the subscript Y denotes the value concerning vortices. The 
following contimity condition must be satisfied fora row of vortices of 
strength T at tiie matching plane g = 


At 6 = 9„ 




-y- = Y . 

Yio yao 


(io) 


- *Y10> 


( 11 ) 


In the similar manner as described for sources/sinks, the constants 

C ,, C ^ can be obtained using conditions (l0) and (ll) as,_ 
ml ' 


Yo N 


1 


■^0 ^ V^20 


Xjji ir , ^\c 


and 


Y ST 

'o 


{_10 _20 _ 

f’ f 
^20 10 

1 _ 


m 






V " , fao fio , V ” ^0 

11““ ^ J 

^10 20 


Using constants Si^^Ujy eqn. (l2) and (l3), we get 


(12) 


(13) 


3$ 


f'' f, 


^ 

Y 38 


/ Y V J_^^oos X (6 . 0 ) d5 - 

" _ J i, -ra f -f m 0 U 


V ‘o f-in ^on 

^ o m=1 0 10 20 


1 “ f ' f, 

^ f J J ^ — 1— COS 

’ "o ^20 ^10 “ 


X ( 9- 9 ) a 5 
m ^ 0 o 


(14) 


and 


'"Y I, « 


=--/y I 

ir ^ ^ , 


1 


’'p* 0 m=1 m 0 10 iji 


sin ( 9 - 8^) d - 


0 

I 


f ' fl 


-Ty^ I 

^ i ° Itt=1 V ^20 ^10 


( 15 ) 
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iVhere the subscript Y denotes the value concerning vortices. The 
following continuity condition must be satisfied for a row of vortices of 
strength r at the matching plane g = 

6= (6;^ 6^) 


\l0 ^20 


(io) 


=^o("t 20 - *T1o’ =r»(8-6o) 


In the similar manner as described for sources/ sinksj the constants 
0 ^ 1 , obtained using conditions (lo) and (ll) as,^ 


Yo N 


-Y. N flJt. 


20 20 


1 -f' f 

f /JO 20 1 X 


"20 "l0 

. o !_! __1__ - :f£j Vfn 

!k!io/ V’ ^0 

10 20 

sing constants 0^^ ,0^2 given by eqn. (l2) and ( 15), v/e get 


r ; _ i r Y T J-fjo.-^ cos X (a/, e ) d 5 

Y o 0 to 20 


J? oa f 

N r y J_ _J0 -J- cos X ( a - e ) d c 

T c m^1 V ^20 - ^ 

_s;/ y sinx (o-ejds^ 

^ C iit=1 V ^0 ^20 10 


( 15 ) 
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Heferriag to fig. 3(a), velocities in 5 and n direction can 
bs \vritten as, 


C = V cos c + w sin e 


C = V sin s + w cos e 
n 


( 16 ) 

(IT) 


With the help of eqns. (l6) and (l7) and using eqns. (7), (8) , (l4) 
and (15), eqns. (2.18) to (2.2l) can be obtained directly. 



iSPPSNDiX - III 


COEEBCTIGII 5DH aiAIiL lEPAETCJRES ®OM "IR EE- VORTEX" MEaI EDW/ 

\'lhea an infinite number of blades is presenb, oce can view the 
flov/ in two vioffs ; 

(i) As being an irrotational, unsteady and non-axisyiiimetric potential 
flow with an infinite number of embedded sheets of vcxrticity shed 
by tte Tariation of lift along the blade span. 

(ii) As a steady, axisymmetric , rotational mean flow (ayeraged in 
©-direction) , with the blade forces represented as a body force 
field. The latter view is taken here. 


The momentum eqxation : 


m " 



(l) 


Where V is the •velocity vector, t is the time variable, p is the 
pressure field, P is the density of the fluid and I is body force per 
unit mass of the fluid, is easily transformed to Crocco's eqve-tion for 
steady flov/ : 

"I X Y = - V H + E ( 2 ) 

where "I is the vorticity •vector api is gi-ven by, 

. ■ J = v^x T ■ 

p 

and H = f + -x- 

where H is the total enthalpy for an isentropio incompressible flow. 
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In (a, p,3) coordinate frame, the a -comp one nt of this eqa. (2) -is, 


(C ?), 


SH 

SR 


(3) 


Since F^, the blade force in a_direction has been assumed to be negligib 
the 000 .( 3 ) can be v»ritten as, 


% ’"p* 


as 




da 


r V 


P* 


-•slGO in the curvilinear coordinate system (a,S,6), we have 


( 4 ) 


p* 


(V X 


9 , 


(— ^)1 
r V ' -1 
p-! 


36 ^ 


1_ 

9a 


(t\7*) H 


''p* 9a 


(5) 


( 6 ) 


as the first term is zero for axisymmetric flow. 


Using eqns. (5) and (6), eqa. ( 4 ) can be written as follows 
3 /■ \ 2 9 4 ** 


v^ 3 ti 


** '■p* 30. 




3H 


r r 


p* 3a 

Srom the Siler’s turbine equation, w© have 

dH = »d (rw^) 


(7) 


( 8 ) 
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Assuming tih tiis flow is a, **siDipl6 SQuilibrium flow** i»6, a flow 
approximately along the a= constant surfaces, u* can be neglected as 
compared to t* or So the eqn.(7) reduces to, 


1 a / n 2 _ _ 1 


2 8a y 


(— ) 




w- 


W* 8 

-liT o wj 


(9) 


^bstituting k= (v^ - into eqn. (9), and neglecting 

2 

the -ceKas of the order of k , as it is assumed that the change in t* from 
V » will be small, we get 


3a 


(•<) = C'*'- ^3 (r w^) 


p* 

a 




or 


h ’’p* 


(10) 


where the subscript h denotes the value at the hub. The volume flow Q 
throu^ the passage of turbomachine can be ejjpressed as follows : 


tip 

Q = / 2^ r V* d Sa 

hub 


(It) 


low using eqn, (10 ) and the definition of k , \ve get 


*^t ® w o 


'h 


Q a' 


r -‘3a 


( 12 ) 


h h::’p* ■ : ; 

where the subscript t denotes the value at the casing and we take 

= 1 ? 



APPfilTDLX IV 


THE COrffiUTER PROGRAJIS 

Two comp'Aer programs were used for the numerical solution of 
the inverse probleia. Ihey are, (i) Por meridional potential duct flov? 
solution, and (ii) 5tor the inverse problem, in which the meridional 
potential duct flow solution is used as an input. 

The computer programs were written in Port ran IV and run on , 
an 7044 computer with approximately 52K words of core stora^. The 
programs are quite general except that sane sections will require changes, 
which will depend on the turbomachine geometry. 

Portran listing of these programs is given in Appendix VI. 

1« Description of Program for Meridional Potential Duct Plow Solution t 

HAIM Program reads and prints various input and output data, solves 
the system of simultaaeoijs linear algebraic equations resulting from the 
finite difference approximation of partial differential eqn. (2.l) and 
calls various sxJbrout ines . 

Description 

Called by MAIM, locates the boundary points, their 
associated intercepts, and ccefficientsG,D,E,P,C 
appearing in eqn. (l) of Appendix I. 

Called by MAIM, for determining positions of a given 
number of streamlines, slope and velocity along them. 
Called by STREAM, for calculating derivative. 


Subroutine Marne 
BMlffTS 

STIiEAt.! 

DBRIV 
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Sx;ibro' .xc ii:»3 llame 

SiOQIH 

iOETRAII jjnpHaiferrt.' 
Pro»3y;3ni. Symbol 
(MArj) 

AIPHA, BEIA 
DX 

D,E,P,G 

ESPS 

EpaiAx 

STSHEP, PISHEP 
IR 

MAI 

jCEtPE 

ITER 

UMAX 


Descript ioa 

C:jlled by SIREAII, for smoothing coordinates of 
streamlines, slope and velocity along thou. 

LISP OPPRIiiaiEAI 7^RD‘iJBLES 
at ioa ; 

Definition 


Grid size of the flovi/ region to be analysed, 
liesh size, ax. 

Coefficients D,E,P,G defined by eqn. (l) in 
Appendix I, and multiplied by C. 

Sum of percentage variations in the functional 
values at grid points from iteration to iteration. 
Converge rue criteria for EpS. 

Shifts in coordinate axes. 

Matrix such that IR(K,J) contains the extreme 
column subscript of a completely interior point 
in row J n.B;ar boiEidary K# 

Matrix such that RIAX (K,J) contains the column 
subscripts I of tie first and last boundary points 
(of types ;I-I7) in row J near boundary 
Matrix giving type of each boundary 
to IV) , shown in fig. l(h). 

Counter on the number of Gauss-Seidal iteration. 
MaxiiQum allowable value of ITER. 



’OTOirraiii Symbol 
"*•''***'’" '■ '' 


Definition. 


I,J 

Grid-point sT±»scripts i,o. 

K 

Subscript in EUX and IH for inner and outer 

boundaries. 

Ii,N 

Number of intervals, in which ALPHA -and BETA are 

divided. 

ITS 

Humber of stream surfaces. 

P3I 

Ilatrix cont.aining values of at each grid point. 

pSEluX 

Maxiinum ralue of 'I' i.e. at the casing. 

IIS’ 

Relaxation factor. 

(Sxibroiitine BNEPTS) 

A,B 

Intercepts near the boundary grid point in BEEA 

and AIPEA directions. 


very small distance. 

(S\abroutine SIRM’H) 

A 

vector containing values of at given number of 

stream surfaces. • 

AAjHB 

vectors containing function and independent variable 

for calculating derivatives. 

AP, BP 

Values of derivatives. 

, MI 

yector oortaiixlng the value of nmi.er of points on 


each stream surface^ 

M,IO 

Counters on the number of points on a stream surface 

SLOPE 

Slope of the meridional streamline. 

TSip 

True Slope of known streamlines. 
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i rp.-grcja Symbol Itefiratioii 

? Meridional velocity along a streamline. 

YRfYZ Velocity components of V. 

XfY Matrix containing coordinates of streamlines. 

(Suoi'outiae SldOO-rK) 

VII,X.l 5 ’iM,S.! Smoothed values of V, X, Y and SlOpE at a point. 

IS'IP Temporary storage for old values used for 

smoothing. 


2, Desorlptlon of Program for the Inverse Problem ; 


ilAII program reads and prints various input data, some output 
data and calls viirious subroutines. 


Subroutine Name 
STRIil 


SSS 


mum 


SiiiiPE 


IIDIUS 


lesGription 

Called by MAII, calculates lengths of various 
streamlines from the inlet edge and non- dimensional- 
ises various input data. 

Punction used for calculating length of a curve 
between any two points. 

Calls a by MidlNjinteipolat es values of specified 
velocity distribution from the ir:put data for 

the required number of points. 

Called by MAIN, calculates initial sh^e of the 

canber 

Function used for selection and/or inteipolation 
for a value of r corresponding to any value of Sg. 
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Stibrouciae laae 


7EEP 


W3MSli 

DUBTYj pmE 
TK/iSHI 

SRHCBN 

7S0URC 

BLiiDE 


Desoription 

Function used for selection and/or interpolation 
for a ralue of corresponding to any value of r. 
Function used for integration between any two 
points. 

Called by MAIU, calculates various variables 
appearing in 1KB solution. 

Called by \iyKBiPE, -used for calculation of derivatives. 
Called by MAIST, calculates induced velocities due 
to vortices. 

Called by MAII, calculates the source distribution 
on stream surfaces. 

Called by MAII, calculates induced velocities due 
to sources and sinks. 

Called by MAII, calculates the camber line shape, 
Blade thickness, estimates v^ from, v ^ and checks 
the convergence of blade profile. 


list OF PRIRCIIAI lAffllABlES 
FOBTRM Iimpiementation 


Program. Symbol 

(maii) 

Al ,11 

APSI 


LefiaLt ion 

Yectors used for different teu^jorary storages . 
Vector contain values of at a number of stream 
suif aces.v 
Blade thickness. 
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pro,':r.-n gyabol 
CLL 


GGI 


caia.io, cETGjc 


CGIQCyGiSfQC 


D 


EiY 


lEDGH 

JJpS 

I?1,If2,]!'1P,E2P,F0 


STL 

ITi® 

irm.i 

!.i 


lffiLA33E 

IS 

IP 

III 

lUIvIBSE 


Definition 

Matrix containing lengths along the camber lines K • 
Specified velocity distribution. 

Induced velocities due to vortices in g and n 
direction. 

Induced velocities due to sources and sinks in C 
and n direct ions. 

Constants from eqn. (4. 1l). 

Dummy rari.able. 
dt/d5 . 

Angle e . 

Functions defined by eqns, (5.8), (5.9) and (5.15). 

Counter' used for nuinber of iterations. 

MaxinMu nimber of terms can be taken for 3.nfinite 
series suamation used for calculating induced velocities 

due to singularities. 

lumber of blades. 

lumber of stream surfaces. 

Number of terms used in Fourier's ApproxL nation. 

Number of points known at each streamline. 

Number of points at which specified velocity 
distribution at each stream surface is given as an 

iDpirt# 
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prOftT':ai S.mbol 

Definition 

P, PP 

Punctions a and a' in eqn. (5.9). 

pSEIAX 

Ilaxiijiuni value of 4* » i.e. at the casing. 

RG 

Coordinate r of the camber line. 

SOURCE 

Source strengths of sources placed on camber line,o , 

3C3EPA 

Length Sg for different points on the camber line. 

SBiPA 

Length of streamline at different points from 


the leading e(^e. 

THEPA 

Coordinate 6 . 

'J}AU 

Punction t . 

VS 

V* 

VST/iR 


VORTEX 

Strength of vortices placed on the camber line, Y , 

V 

Meridional potential flow velocity. 

X 

Z-coordinate of streamlines. 

XX 

value of 5/Jl for the points at vdiich velocity 


distribution is specified as an input. 

Y 

Coordinate r of streamlines. 


( Sti) routine STRIM) 

mis ^tip* 

RTIPS * 

(Stib^routijtiG SHi®B) 

Oireulatiott, r defined by eqn. (2.42). 

Coordinate r of streamlines. 

Slope of the camber line. 


TMapS 
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Pi-o;'a"-an Symbol 

( Subrouti ne 'wTCBiiT R) 

ai 

G1,32 

t, 

(Subroufciiie DSRIV) 
YS 

(SubroTjfeine 

/jaiG 

/gaic 

BH 

CON 

FIS’ 

BOY® 

®OYFE 

VjOu 

(SxJbrmtiae SRBEBN) 
ilQG . 

iiJQO 


Befinitioa 


Function defined by eqn. . (2.12). 

Functions defined by eqn. (2.1 1). 

Function defined by eqn. (2. 14). 

Fu.iction defined by eqn. (2. 10 ) and its first and 
second derivatives. 


Derivative. 


defined by eqn. (2.20 ). 

, defined by eqn. (2.2l). 
defined by eqn. (5.6). 

Convergence criteria for infinite series 


simmation. 

fft f! Wff f p ),• defined by eqns. (2. 20) 
*^^10 20-^^^^ 10 20 

or (2.2T). 

f^f ’ , defined by eqn. (5 .14) . 
f'/f'» defined by eqn. (5*12). 

iibsolute value of last term of the series to be 
added in the series sunmat ion. 


I . defia.6<i05y eqrx* (2.18). 

. qc ^ ■ 

j , defined by eqn. (2.l9). 
qc’ 


iJQGYX 



Sysisol 


Itefinitioa 


BK 

defined by eqn.(5.6). 

/I T 

V'kT*** '*' “*■* 

Indiced velocity due to cascade of sources, defined 

E3 

by eqn. (2.55). 

E, define by eqn. (2.37). 

m 

percentage change in value of source strength at 

EREOR 

some point. 

Cumulative percentage change in value of source 

POYF 

strength at all points on a camber line, 
f^f, defined by eqn. (5.13) . 

EPOYF 

fVfj defined by eqn. (5.10 ). 

ITER 

Counter for iteration in the solution of integral 

T 

eqn. (2.40). 

Function defined by eqn. (2.26). 

(Subroutine ■’JSOURC) 

Symbol used aie saine as already defined above. 

(SubroiAine BIiADE) 

m 

K, defined by eqn, (2.49). 

Msm 

a, defined by eqn. (2.52). 

EREOB 

cumulative percentage variation at aH poitxfis on 

m 

a camber line. 

Function inside the integral in eqn. (2.49)^. 

MH 

^ , defined by eqn. (2.5l). 

&iaa \ 

Circulation, T defined by eqn. (2.42). 


product (r t%)« 

SAWH 

■ Blaae tMolEBS® of tHe previous iteratioo. 

STHESA 

e ^oordiaat « “f P^'^ous iteratloa. 



ffii aiDix V 


3A”A fob the mBRICiXL EXJffiLE 


f jsrious i:touc data used for the numerical exaa^ie given 


heruj are given below 

. iJiafcs of these variables axe same as used 

in Bortri-m progra;:i. 


Meridional lotential. 

duct flow Solution : 


Stream function values at various stream surfaces are: 


0.0(1), 0.5345(11), 1.069(111), 1.6035(1?), 2.138(7). 

AEPlie’i 

80 ram. 

31'® A 

100 mm. 

ijpaiiAx 

0.01 

ITMAX 

400 

M 

40 

H 

50 

IS 

5 

RP 

Rf,1.8 

Inverse iroblem 


APSI 

Same as A given above. 

anh 

Total oamb® line length In aot^dinimsiotea fom are 
1.75(1), 1.33(11), 1.05(:il), 0.855(17), 0.65(/). 

:m A 

Maodmun number of terms for infinite series sunmation: 

IffilADB 

1,8 

m 

17 

lUMBm 

13 

PSIMAX : : r 

2.138 (nort-dimensional) 


r^j^p,6Tmm. 



r 


* * ®1 s 


® i «■ 
, ® « 
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ft » , ■ i ■ ' V ^ > -Al ® 

!». S » I© t< ® « * ■« » ■«■ ® m ^ ^ m ■9 0 9 ® 9 9 ® 9 


4 ix it >lc ix 'X ❖ 'X ❖ 'X t s^c :^c ❖ V 5§5 • 


! I :.J F , ii 1) ! ( 2, 5 1 ) t I T YPF ( 26 , 36 ) , D ( 2 6 , 3 6 ) , ■:. ( 26 , 36 ) , 


1 I i :.l ' 2 j *> '’> ' » -' { Afi » .3 1> ) j P S i ( 4 i t 5 L ) 

' ’-.r'V’:' I ''*f ■«»:■■;'? IT-1A.X ,ALPHAtBrTA,E PSMAX,RF 
;■' A I : F 1 r , A, IIP AX , ALPHA, B-TA, LPSHAK ,RF 
.Pl-^'if.l 
: ; f* 1 ■= •'■. ■^- 1 
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■ 1 P'-' .P, ' > 
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■:p' I P 
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* 1. J m ''■ ) 


■it,* 1 

I = !,.Mp.i 


^ \:t i 
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1 

P S i { 
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OC f 

I=2Vf'Al , 

2 ■ 

' 'PS iC 

•I ', 1 ) '” ' .* • '- > ■ 


1.V u ' ..f 

J= -J.fi,"*,)! 

3 

' . P^$I ( 

13, J 5:=PSPiAX 


OC 4 

J = 2.6,61 ..- 

4 

. PSl I 



' P' 3 '*■■'■' P-v 
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.6 ' 
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. -OC '■■■■■ PA.'- ■ .' „ 
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)) A^iC C?TEPM^I NATI ON OF ITS TYPt, 
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A 1 i. 1 
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r-’SI VaLs. SP':CIFI£D uM I-ACH STRtA'’! SURFACE. 


{ I ; . A 


A! 'A? =A 

C' 

::.,.X/4 

J 

% 

-•»(2) 
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r s| j Pl'U' i V ( ^ J 

V7=-{ U7P7P«AP 


:.yXi</0'X,: 


A'7r-2 i*OX 


A A'F 3' ) ®'i2 


( 3^5siA^.{ IX? J+i ):*X l,-D 
IF{ IMAXX 1, 772) -GO . I , j ) IGO TO 'fil, 
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